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Abstract 

A first measurement of the cross section of the process e + e~ — > Z77 is re- 
ported using a total integrated luminosity of 231 pb" 1 collected with the L3 detec- 
tor at centre-of-mass energies of 182.7 GeV and 188.7 GeV. By selecting hadronic 
events with two isolated photons the e + e~ — > Z77 cross section is measured to 
be 0.49lg;?? ± 0.04 pb at 182.7 GeV and 0.47 ± 0.10 ± 0.04 pb at 188.7 GeV. The 
measurements are consistent with Standard Model expectations. Limits on Quartic 
Gauge Boson Couplings a /A 2 and a c /A 2 of -0.009 GeV" 2 < a /A 2 < 0.008 GeV" 2 
and -0.007 GeV" 2 < a c /A 2 < 0.013 GeV" 2 are derived at 95% confidence level. 
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1 Introduction 



The LEP centre-of-mass energy (y/s) for e + e~ collisions has now exceeded the W pair and Z 
pair production thresholds, allowing the study of triple gauge boson production processes such 
as e + e~ — > Z77 and e + e~ — > W + W~7. Measurements of these processes give a new insight into 
the Standard Model of electroweak interactions (SM) 1J. The possibility of these triple gauge 
boson production processes proceeding via s-channel exchange of a fourth boson provides a 
probe of quartic gauge boson couplings (QGC). Such measurements were recently performed 
for the e + e~ — ► W + W~7 process j2]. 

This letter describes the first measurement of the cross section of the process e + e~ — > Z77 
followed by the hadronic decay of the Z. In the SM this process occurs by radiation of the 
photons from the incoming electron and/or positron, corresponding to a total of six diagrams, 
three of which are presented in Figure No QGC contribution is predicted at the tree level 
and Z77 events are sensitive to anomalous QGC jSHS] contributions, as shown in Figure ^ 

The measurement uses data collected with the L3 detector jH] at LEP in 1997 and 1998 
at average centre-of-mass energies of y/s = 182.7 GeV and y/s = 188.7 GeV corresponding to 
integrated luminosities of 55 pb -1 and 176 pb -1 , respectively. These energies are respectively 
denoted as 183 GeV and 189 GeV hereafter. 

The e + e _ — ► Z77 — > qqyy signal is defined by three phase space cuts: photon energies above 
5 GeV, photon angles with respect to the beam axis between 14° and 166° and invariant mass 
of the primary produced quarks within a ±2Pz window around the Z mass, Tz being the Z 
width. The KK2f Monte Carlo (MC) program [2] predicts signal cross sections of about 0.4 pb 
at both energies. 

2 Event Selection 

The selection of events satisfying the signal definition given above is optimised using hadronic 
events generated at y/s = 189 GeV with the KK2f MC program and at y/s = 183 GeV with 
the PYTHIA5.72 |S] MC program. Events from these MC programs failing the signal def- 
inition are termed QCD77 background. Other background processes are generated both at 
y/s = 183 GeV and = 189 GeV with the MC programs PYTHIA (e+e - Ze+e" and 
e+e" -»■ ZZ), KORALZ4.02 P (e+e" -> t + t-( 7 )), PHOJET 1.05c TDJ (e+e" -»■ e+e'qq) and 
KORALW 1.21 ^T] for W + W~ production except for the ei> e qq' final states which are generated 
with EXCALIBUR [T2j. Additional background sources are found to be negligible. 

The L3 detector response is simulated using the GEANT 3.15 program |13| . which takes 
into account the effects of energy loss, multiple scattering and showering in the detector. Time 
dependent detector inefficiencies, as measured in each data taking period, are reproduced in 
these simulations. 

The selection of e + e~ — > Z77 — > qq77 candidates from balanced hadronic events with two 
photons and little energy deposition at low polar angles is based on photon energies and angles 
together with the invariant mass of the hadronic system. The photon energy and angle criteria 
follow directly from the signal definition whereas the invariant mass of the hadronic system 
is required to be between 74 GeV and 116 GeV. The main background after these selection 
requirements is due to the radiation of two initial state photons with a hadronic system failing 
the Z signal definition criteria. The boost (3% of the recoiling system to the photons, assuming 
its mass to be the nominal Z mass, is on average larger for these background events. Candidate 
events are hence required to have (3z < 0.64 at y/s = 183 GeV and j3z < 0.66 at y/s = 189 GeV. 
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Another class of background events is the so called radiative return to the Z, where a photon 
in the initial state is emitted, bringing the effective y/s to the Z resonance. A Z boson is then 
produced decaying into a hadronic system with an electromagnetic energy deposition (photon, 
misidentified electron or unresolved 7r°) faking the least energetic photon of the signal selection. 
These events are rejected by an upper bound on the energy Ej of the most energetic photon and 
a lower bound on the angle u between the least energetic photon and its closest jet. Numerically 
Ej < 67.6 GeV at ^ = 183 GeV and E\ < 70.7 GeV at = 189 GeV with u > 17° at both 
the energies. 

Data and MC distributions of the selection variables are presented in Figure EJ for y/s = 
189 GeV, where selection criteria on all the other variables are applied. A good agreement 
between data and MC is observed. 

3 Results and Systematic Uncertainties 

The application of the selection procedure described above yields the signal efficiencies and 
selected data and MC events summarised in Table 1. A clear signal structure is observed in the 
recoil mass spectra of the two photons, presented in FigureElfor the two centre-of-mass energies 
under study. The e + e~ — > Z77 — > qqyy cross section at y/s = 183 GeV and y/s = 189 GeV is 
then determined from the number of events selected and the efficiency and background estimates 
from MC. 



(GeV) 


183 


189 
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0.49 


0.51 


Data 


12 


36 


MC 


13.4 


39.2 


Z77 


10.6 


32.6 


QCD77 


2.7 


6.0 


Other 


0.1 


0.6 



Table 1: Yields of the e + e~ — > Z77 — > qqyy selection at y/s = 183 GeV and 
y/s =189 GeV. The signal efficiencies e and the number of expected events for data 
and MC are given. 

Systematic uncertainties on the cross section measurements are listed in Table 2. They 
include uncertainties arising from the signal MC statistical error of 3% at y/s = 183 GeV and 
5% at y/s = 189 GeV. The uncertainty on the accepted background due to MC statistics are 
7% and 17% at y/s = 183 GeV and y/s = 189 GeV, respectively. The calorimeter energy scale 
uncertainty is estimated by varying electromagnetic energies by ±1% and hadronic energies by 
±2%. Selection procedure systematics are obtained from the effect of removing each of the selec- 
tion criteria. From a comparison of the KK2f and PYTH1A cross sections for hadronic events 
with a hard photon, a ±15% uncertainty on the background normalisation is conservatively 
estimated and propagated to the measured cross section. Uncertainties on signal efficiencies 
are estimated by comparing KK2f with PYTHIA and GRACE MC program predictions 
and are found to be negligible. The systematic error assigned to the reweighting procedure at 
y/s = 183 GeV is estimated by applying the procedure to a sample generated with PYTHIA at 
y/s = 189 GeV and comparing the corresponding cross section result to the previous one. 

The cross section results are: 
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Table 2: Systematic uncertainties Aa on the e + e — > Z77 — > qq77 cross section at 
v/i = 183 GeV and = 189 GeV. 



a e +e-_z 77 ^qq 77 (183 GeV) = 0.34±g;J| ± 0.03 pb (SM 0.396 ± 0.005 pb) 
cr e + e -_ +Z77 ^ q q 77 (189 GeV) = 0.33 ± 0.07 ± 0.03 pb (SM 0.365 ± 0.003 pb), 

where the first uncertainties are statistical and the second systematic. The values in paren- 
theses denote the SM expectations calculated from the KK2f MC with its default set of input 
parameters. The error on the predictions is the quadratic sum of the MC statistical error and 
the theory uncertainty estimated as suggested in Reference [I] . The measurements are in good 
agreement with these predictions. These results are also presented in Figure 0] together with 
the expected evolution with -y/s of the SM cross section. 

Scaling the measured cross sections for the Z hadronic branching ratio gives the e + e~ — > Z77 
cross section at y/s = 183 GeV and ^fs = 189 GeV: 

<7 e + e -^z 77 (183 GeV) = 0.4918;?? ± 0.04 pb 
a e+e -^ Z77 (189 GeV) = 0.47 ± 0.10 ± 0.04 pb. 



4 Limits on Quartic Gauge Boson Couplings 

Anomalous QGC contributions to Z77 production via the s-channel exchange of a Z are de- 
scribed by two additional terms of dimension six in the Lagrangian [3]: 

£\ = -^a c F w F^W?-W a , 

where a is the electromagnetic coupling, is the field strength tensor of the photon and 
Wfj, is the weak boson field. For Z77 the third component of W^, Z^j cos^iy, is relevant. 
The parameters ao and a c describe the strength of the QGC and A represents the scale of 
the New Physics responsible for the coupling. £g and £{§ are separately C and P conserving 
and no CP violating operators contribute to the anomalous ZZ77 vertex. A more detailed 
description of QGC has recently appeared [5 . While indirect limits on the QGC were derived 
from precision measurements at the Z pole ^H], studies of Z77 and W + W~7 production probe 
the quantities ao/A 2 and a c /A 2 in a direct way. The e + e~ — > Z77 process is expected to have 
higher sensitivities than e + e~ — > W + W~7. This is due to an extra factor of 1/ cos 4 ^ in the 
QGC cross section, to the larger SM cross section and data statistics and to the smaller number 
of SM diagrams jlj. 
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QGC are expected to manifest themselves via deviations in the total e + e" — > Z77 cross 
section, as presented in Figure 0] As the Z77 production occurs in the SM via t-channel 
diagrams, the three body phase space favoured in the QGC mediated production is different, 
in particular resulting in a harder spectrum of the least energetic photon [5.. Figures and 
EJd compare these reconstructed spectra with the predictions from signal and background MC 
at y/s = 183 GeV and y/s = 189 GeV. The expectations for an anomalous value of ao/A 2 or 
a c /A 2 are also shown. These QGC predictions are obtained by reweighting each SM signal MC 
event with the ratio W(0, ao/A 2 , a c /A 2 ), a function of its phase space Q derived from the two 
photons and the Z mass and the values of the couplings: 

vwfo /A2 /a^ \M SM {n) + M QGC {n, a /A 2 , a c /A 2 ) | 2 
VV{\l,a /A ,a c A ) — -— . 

M.sm denotes the SM matrix element and M.qgc the QGC one, both calculated analytically [I]. 
Possible extra initial state photons are taken into account in the calculation of Q. 

A simultaneous fit to the two energy spectra is performed leaving one of the two QGC free 
at a time, fixing the other to zero. The SM predictions in the fit procedure are reweighted as 
described above, yielding the 68% confidence level (CL) measurements: 

a /A 2 = 0.001 ± 0.004 GeV" 2 and a c /A 2 = 0.003 ± 0.005 GeV" 2 , 

in agreement with the expected SM value of zero. A simultaneous fit to both the parameters 
yields the 95% CL limits: 

-0.009 GeV" 2 < a /A 2 < 0.008 GeV" 2 and - 0.007 GeV" 2 < a c /A 2 < 0.013 GeV" 2 , 

as shown in Figure El A correlation of —35% is observed. The experimental systematic un- 
certainties and those on the SM e + e~ — > Z77 — > qqyy cross section predictions are taken into 
account. 
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Figure 1: Three of the six SM diagrams contributing to e + e~ — > Z77 production. 
The other three SM diagrams are obtained by crossing the photon lines. A possible 
anomalous QGC diagram is also shown. 
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Figure 2: Distributions of (a) the invariant mass M q q of the hadronic system, (b) 
the boost /?z of the reconstructed Z boson, (c) the energy Ej of the most energetic 
photon and (d) the angle u between the least energetic photon and the nearest jet. 
Data, Z77 and background MC are displayed for ^/s = 189 GeV. The arrows show 
the position of the final selection requirements. All the other selection criteria are 
applied for each plot. 
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Figure 4: Evolution of the e + e~ — > Z77 — > qq77 cross section with the centre-of- 
mass energy. Signal definition cuts described in the text are applied. The width of 
the band corresponds to the error that arises from MC statistics and theory uncer- 
tainty, estimated to be 1.5%. Dashed and dotted lines represent QGC predictions. 
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Figure 5: Energy of the least energetic photon for (a) y/s = 183 GeV and (b) 
y/s = 189 GeV. Data, Z77 and background MC are displayed together with QGC 
predictions. 
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Figure 6: Two dimensional contours for the QGC parameters Oq/A 2 and a c /A 2 
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